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CLOUD COMPUTING AND HYPERBOLIC VORONOI DIAGRAMS ON 

THE SPHERE 

PAVEL BLEHER AND CAROLINE SHOURABOURA 

Abstract. In this work we study the minimization problem for the total distance in a 
cloud computing network on the sphere. We prove that a solution to this problem is given 

f— ^ in terms of a hyperbolic Voronoi diagram on the sphere. We present results of computer 

P\J simulations illustrating the solution. 

(^ 1. Introduction 

We consider a network 
^ S = {Si, . . . , Sn} 

I of computer data centers on the earth surface. We assume that the earth surface is rep- 

^ resented by a sphere A/j of a radius R > 0, and each data center Sm is described by its 

geographical coordinates, the latitude (p^ and the longitude 9^- The data centers provide a 
variety of cloud computing services to a big set of users 



cd 



^ U = {Ui,...,Ur,}. 

,__! Our main goal will be to solve the minimization problem for the total distance T) between the 

lO data centers and the users. We assume that each user Uk is assigned to a data center Sji^k), 

^ and the total distance is the sum of the distances between the users and the corresponding 

• data centers, 

cn 

O "^ 

CM V{j) = Y,d{S,^k),Uk), (1.1) 

• ^ where j = j{k) is the function assigning the user Uk to the data center Sj(^k)- 

^ For any two points x, y on the sphere we denote d{x, y) the distance between x and y 

on the sphere, that is the length of the arc of the great circle connecting x to y. We will 
consider the problem of the minimization of the total distance under certain constraints. 
We will assume that each data center Sj has a capacity Cj which characterizes how many 
users it can service. The assumption of the limited capacity for the data centers leads to the 
minimization problem with constraints: Find an assignment function jo{k) such that 

min V{j)=V{jo), (1.2) 

where 

J(C) = {j=jW|'V.n = #{*: jW = "}<C„. for m=l,...,n}. (1.3) 
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Observe that Nm = i^{k : j{k) = m} is just the number of users assigned to the data center 

Since the number of users is usually large, it can be useful to consider a (nonnegative) 
measure /i of users on the sphere. In this case, the total distance functional looks like 

Vij)= [ diS,^u),U)dfx{U), (1.4) 

where j{U) is the assignment function of the users to the data centers. Respectively, con- 
straint ( 1.3[ ) takes the form 

J{C) = {j=j{U)\fi{U: j{U) = m}<C^ for m = l,...,ri}. (1.5) 

The main result of this work is Theorem 2.1, in which we prove that the solution to 



minimization problem (1.2) is given by a hyperbolic Voronoi diagram on the sphere. The 
classical Voronoi diagrams have numerous applications to science and technology (for many 
of these applications see the review paper [Ij by Aurenhammer, the lectures |6j by M0ller, 
the collection of papers |1], the monograph [8] of Okabe, Boots, Sugihara, and Chiu, and 
references therein) . For applications of Voronoi diagrams to cloud computing see the recent 
paper [U] of Shouraboura and Bleher and references therein. In our Theorem 2.1 we show 
that the presence of the constraints leads to an extension of the classical Voronoi diagrams 
to hyperbolic Voronoi diagrams on the sphere. It is interesting to notice that in a completely 
different setting similar diagrams on the plane or in the 3D space appear in crystallography 
as the Johnson-Mehl model of the crystal growth [5]. We formulate and prove Theorem 2.1 
in Section |2j In the subsequent section we illustrate Theorem 2.1 by results of computer 
simulations. 

2. Solution of the constrained minimization problem and hyperbolic 

Voronoi diagrams on the sphere 



In this section we will solve minimization problem (1.2). We will consider a general 
measure fi of users. For a given measurable assignment function j{U) on A/j, we introduce 
the sets 

amij) = {U eAr: jiU)=m}, m = l,...,n, (2.1) 

so that <7m{j) is the set of users assigned to the data center Sm- The sets {o-mij)} are disjoint 
and they form a measurable partition tt{j) of the sphere Ar, so that 

n 

n{j) ■■ Ar=[_\ amU)- (2.2) 

m=l 

The functional Vi^j) reduces to 

n „ 

^0') = E/ d{Sm,U)dfi{U), (2.3) 

and constraint (1.5[) to 



J{C) = {3=3{U)\fi{aUj))<Cm for m = l,...,n}. (2.4) 

Let us begin with the unconstrained minimization problem. In the absence of the capacity 



constraint, minimization problem (1.2) can be solved by assigning each user U to its closest 
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data center Sj, so that Sj^^u) is the closest data center to U. Geometrically, we partition the 
sphere A^j into the cells {am} of the Voronoi diagram V on the sphere with the points {Sm}- 
The cell am is defined as the set of points on Ar which are closer (or at the same distance) 
to Sm than to other S'/'s, that is 

Cm = {x E Aji : d{x, Sm) < d{x, Si) for all / 7^ m}. (2.5) 

Observe that the cells {am} are convex spherical polygons on the sphere A^^. With the 
Voronoi diagram we associate a graph Fy. The vertices of the graph Fy are the vertices of 
the polygons {am}, and the edges of the graph Fy are the sides of the polygons {am}- The 
Delaunay triangulation, associated with the Voronoi diagram, is the dual graph Td to Ty, 
with vertices {Sm} and edges connecting vertices 5*^, Si if and only if the cells am, cri have 
a common side. 

Thus, in the absence of the constraint, we assign to each data center Sm all users in the 
cell am of the Voronoi diagram. To describe the assignment in the presence of the constraint 
we introduce hyperbolic Voronoi diagrams on the sphere. 

Suppose that to each point 5*^ a number dm is assigned. Then the hyperbolic Voronoi 
diagram 

V{d), d= {di,...,dn), 
with the parameters {dm} and the points {Sm}, is defined as follows. The cell am of the 
point Sm is defined as 

o'm = {x & A/j : d{x, Sm) + dm < d{x, Si) + di for all / 7^ m} . (2.6) 

Observe that on the curve 'jmi separating two neighboring cells am, cri we have the equation, 

d{x,Sm) + dm = d{x,Si) + di, (2.7) 

hence jmi is a part of the spherical hyperbola on the sphere Apt- A vertex v of the hyperbolic 
Voronoi diagram is a point which belong to three or more cells. The graph Tv{d) of the 
hyperbolic Voronoi diagram V{d) consists of the vertices {v} and the edges {e}, which are 
the curves {■jmi} separating neighboring cells. 

A solution to the constrained minimization problem can be obtained as follows. We will 
take a natural assumption that the total capacity of all data centers is not less than the 
number of the users, 



2_^c„ 



m=l 



> Mi?- (2.8) 



Otherwise, it would be impossible to service all the users. 

Theorem 2.1. For any measure fi, a minimizer Jq{U) to constrained minimizing problem 
(1.2) exists, which can be obtained as follows. There exist numbers {di, . . . ,dn) such that 
the minimizer jo is obtained by assigning all users in the cell am of the hyperbolic Voronoi 
diagram V{d) to the data center Sm, so that 

Jo{U) = m if and only if U E am- (2.9) 

Proof. For any e > 0, we will introduce a regularization fis{U) of the measure ^{U), and we 
will first prove Theorem 2.1 for the measure fie{U) and then we will take the limit e — )• 0. 
The measure fie{U) is defined as follows. Consider the function 

^,(x) = e-^"/^ (2.10) 



on the real line, and for a given point Uq G Ar, introduce the probability density fe{U; Uq) 
on the sphere Ar as 

fe{U- f/o) = ^9e {d{U, f/o)) > (2.11) 

where 



ZiE)= / g,{d{U,Uo))dU, (2.12) 

Then, as e — )■ 0, the function /^(f/; f/o) converges to a 5-function at the point Uq. Consider 
the function 



p,{U)= / MU;Uo)dfi{Uo)>0, (2.13) 

JAu 

and define the measure HeiU) as 

dfi,{U) = Pe{U)dU. (2.14) 

Then in the weak sense, 

lim//e = /i, (2.15) 

so that for any continuous function f{U), 

hm / f{U)dfi,{U)= [ f{U)dfi. (2.16) 



Consider the minimization problem 

min V,{j)=V,{j,), (2.17) 

for the functional 

V,{j)= [ d{Sj^u),U)dfie{U). (2.18) 



Lemma 2.2. yl minirnizing function j^ for minimization problem (2.17) exists. 

Proof. Let 

inf P,(j)= limP,(j,^). (2.19) 

Consider the partitions vr(j^). By the compactness argument, there exists a weakly converg- 
ing subsequence 

lim 7r(j^) = 7r(j,). (2.20) 

r— >oo 

Since d{U,Uo) is a continuous function of (f/, f/o) and PeiJJ) is a continuous function of f/, 
this implies that j^ solves (2.17). D 



Lemma 2.3. Let j^ he a solution to (2.17) and 



m=l 

the corresponding partition of A^. Then for any m ^ I, 



Ue) ■■ ^R=U ^-(•^■-) (2.21) 



essinf [d{U, Si) - d{U, S^)] > ess sup [d{U, Si) - d{U, Sm)]. (2.22) 



Proof. Suppose, for the sake of contradiction, that (2.22) is false. Then there exist sets 
A C (Jiije), B C (Jmije) and numbers a > b such that 

(1) d{U, Si) - d{U, S^) > a for all U e A, 

(2) dlu, Si) - dlu, Sj) < b for all U eB, 

(3) fie{A) = fie{B) > 0. 

Let us reassign all users U E A from Si to Sm and all users U E B from 5*^ to 5*;. This 
reassignment, j^, will not violate the constraints and it will decrease the total distance 
functional V^j). Indeed, 



^.a)-^.(i 



[d{U, Si) - d{U, SJ] dfi,{U) - / [d{U, Si) - d{U, SJ] d^,{U) 



B 



> aii,{A) - biie{B) = (a - b)ii,{A) > 0, 



hence ^(je) is not the minimum of T>{j). This contradiction proves (2.22). 



(2.23) 
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It follows from (2.22) that there exist sets {Em, m = 1, . . . ,n} such that 

(1) Err, C amije) ioi 171 = 1, . . . , U, 



(2) /i£ (amiJe) \ Em) = for m = 1, . . . , n. 



(3) for any m ^ I, 



inf [rf(f/. Si) - diU, Sm)] > sup [d{U, Si) - d{U, Sm). 



Consider the closure Fm of each set En 



Fm = Em, m = l,...,n. 



Then, since {cmije)} is a partition of the sphere An, we have that 

n 
Ai? = U Em. 



(2.24) 
(2.25) 
(2.26) 



rra=l 



By continuity, (2.24) holds if we replace Em, Ei for Em, Ei, respectively, that is for any m ^ I, 

(2.27) 



inf [d{U, Si) - d{U, Sm)] > sup [d{U, Si) - d{U, Sm). 
Moreover, if the intersection of Em and Ei is nonempty, 

-fml = EmUEiJ^iD, 

then 

inf [d{U, Si) - d{U, Sm)] = sup [d{U, Si) - d{U, Sm) = dmi. 
ueFm ueFi 



This implies that for U G 'jmi. 



(2.28) 
(2.29) 

(2.30) 

Thus, 7m/ is a part of hyperbola (2.30) on the sphere Ar. Since the function d{U, Sm) is an 
analytic function of U outside of the point 5"^, it follows that 'jmi consists of a finite number 

(2.31) 



d{U, Si) - d{U, S„ 



d 



ml. 



of arcs of hyperbola (2.30). Moreover, the boundary of Em is equal to 

dEm= [jiEmnEi), 

Ij^m 

hence dEm consists of a finite number of arcs of hyperbolas on the sphere A^^. 



Assume that for some m the sets 'jmi and 'jmp intersect at some point v. Then v belongs 
to the boundaries of Fi and Fp, hence 7/p passes through v, that is f is a triple point (or 
higher). By taking U = v in (2.30), we obtain that 

(^ml ~r dip ~r dtpm U. 

Define the one-chain 

Cl = / ^dmllml- 
ml 



(2.32) 
(2.33) 



Then by ( l2.32[ ), 

dCi = 0. (2.34) 

Since the first homology group of the sphere Aji is trivial, this implies that there exists a 
zero-chain Co such that 

dCo = Ci. (2.35) 

In other words, there exist numbers {dm} such that 

dmi = dm — di. (2.36) 



By (2.30), this proves (2.7) and hence Theorem 2.1 for the distribution /i^ of users. Let us 
consider the distribution /i. 

The numbers {dm} are defined uniquely up to a common additive constant. To fix the 
constant, let us put di = 0. Then the other numbers, dm = dm{£), m > 1, are bounded by a 
constant independent of e. Let us take a converging subsequence 

(i™, = lim c?™,(e„), p = l,2, ...; m = l,...,n. (2.37) 



lim dmisp), 



This defines the numbers {dm}- Consider the corresponding hyperbolic Voronoi diagram 
V{d). Since V{d{e)) gives a solution to the minimization problem for the users distribution 
fie, we obtain that V{d) gives a solution for the distribution /i. This proves Theorem 2.1 D 



Remark. Observe that, in general, the hyperbolic Voronoi diagram solving the minimiza- 
tion problem is not unique. Consider, for instance, a discrete distribution of users 



N 



dii{U) = ^6{U -Uk)dU. 



(2.38) 



k=l 



Then, if there are no users on the edges of the Voronoi diagram, then we can shift the 
numbers di, . . . ,dn a little, without changing the number of users in the Voronoi cells. On 
the other hand, the solution to the minimization problem is unique for any users distribution 
which has an absolutely continuous component with a positive density, like /i^. This explains 
why we used the distribution /i^ in the above proof, establishing the theorem first for fj,^ and 
then passing to the limit e — t- 0. 

3. Computer simulations of the minimizer 
Consider a minimizer jo of ( |l-2 ). We say that constraint (1.5) is active on a server Sm if 

fi{U : jo{U) =m} = Cm, (3.1) 



and it is inactive on Sm if 



fx{U : joiU) = m}<C„ 
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(3.2) 



If all constraints are inactive then the minimizer jo is the usual Voronoi diagram on the 
sphere. 

In the computer simulations we begin with the usual Voronoi diagram. To that end, we 
construct the convex hull "H of the set S. This can be done in 0{n\nn) operations. The 
convex hull gives the Delaunay triangulation, and the spherical bisectors of the edges of the 
Delaunay triangulation give the edges of the Voronoi diagram (see, e.g., [2], [3], [Z])- We 
check the condition 



Ai(o-m) < C„ 



(3.3) 



for the Voronoi cells am, rn = 1, . . . ,n. 
minimizer is the usual Voronoi diagram. 



If the condition is satisfied for all m's, then the 




Figure 1. The solution to the constrained minimization problem with 10 
data centers, Seattle, Atlanta, New York, Phoenix, San Francisco, Denver, 
Houston, Chicago, Boston, and Miami. The Voronoi Diagram on the sphere 
with these 10 points is shown by the cyan line. We impose two constraints 
shown in Red: in San Francisco the capacity is equal to 15 million users and 
in Atlanta it is 20 million. 



In application to cloud computing, it is plausible to assume that most of the servers have a 



big capacity, and condition (3.3) violates only on a set Sq of isolated S'm's. To set a notation, 
let 

So = {Sm, m G Jo}, (3.4) 

where /q C {1, . . . , n}. In this case we are looking for parameters {dm, rn G Jq} such that 



^{Om{d)) = Cm, m G Jo, 



(3.5) 



where (Tm{d) is the cell of the hyperbolic Voronoi diagram V{d), with given dm when m ^ I[ 







and dm = when m ^ Iq. We solve equation (3.5) numerically for each m G Iq. We 



assume that the capacity of the neighboring vertices to Sm ^ <So is big so that the resulting 



hyperbolic Voronoi diagram V{d) satisfies condition (3.3) for all m. 

In Fig. [T] we illustrate the numerical solution by an example. In this example, we con- 
sider 10 data centers located in Seattle, San Francisco, Phoenix, Denver, Chicago, Houston, 
Atlanta, Boston, New York, and Miami. We assume that the distribution of users is propor- 
tional to the population distribution. The blue lines show the spherical Voronoi diagram. We 
further assume that the data centers in San Francisco and Atlanta have a limited capacity 
equivalent to the population of 15 and 20 million, respectively. We calculate the hyperbolic 
Voronoi diagram under these constraints. The change from the usual Voronoi diagram to 
the hyperbolic one is shown in Fig. [T] in red. 

4. Conclusion 

In this work we studied the minimization problem for the total communication distance 
in a computer cloud under the condition of restricted capacity of the data centers. Our 



main result is Theorem 2.1, which shows that a solution to the minimization problem is 
given by a hyperbolic Voronoi diagram constructed on the data centers Si, . . . , Sn- The 
parameters di, . . . ,dn of the hyperbolic Voronoi diagram can be found from the condition 
that the number of users in each cell am of the diagram does not exceed the capacity of the 
corresponding data center Sm- 

Although we discuss the apphcation to the computer cloud only, the hyperbolic Voronoi 
diagrams on the sphere can find other important applications. We can mention the problem 
of location of air-bases [8], the distribution of facilities in global Internet companies like 
Amazon.com, the distribution of the telecommunication centers for mobile telephones in 



global telephone companies, data collection centers, and others. Also, theorem 2.1 can be 



extended, under appropriate conditions, to Riemannian manifolds in dimension 2 and higher. 
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